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Gene expression is controlled by the joint effect of (i) the global physiological state of the cell, in 
particular the activity of the gene expression machinery, and (ii) DNA-binding transcription factors 
and other specific regulators. We present a model-based approach to distinguish between these two 
effects using time-resolved measurements of promoter activities. We demonstrate the strength of the 
approach by analyzing a circuit involved in the regulation of carbon metabolism in E. coli. Our 
results show that the transcriptional response of the network is controlled by the physiological state 
of the cell and the signaling metabolite cyclic AMP (cAMP) . The absence of a strong regulatory effect 
of transcription factors suggests that they are not the main coordinators of gene expression changes 
during growth transitions, but rather that they complement the effect of global physiological control 
mechanisms. This change of perspective has important consequences for the interpretation of 
transcriptome data and the design of biological networks in biotechnology and synthetic biology. 
Molecular Systems Biology 9: 634; published online 22 January 2013; doi:10.1038/msb.2012.70 
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Introduction 

Bacterial cells continuously adjust gene expression in response 
to challenges from their environment. These adjustments 
involve transcription factors that sense metabohc signals and 
specifically activate or inhibit target genes. Several hundreds 
of transcription factors have been identified in E. coli (Keseler 
et aU 2011) : while some respond to a particular stress and have 
only a few targets, others coordinate the expression of 
hundreds of genes across a variety of cellular functions. 
Well-known examples of the latter are global regulators of 
transcription, such as Crp, Fis, and RpoS (a^) (Gottesman, 
1984; Martinez-Antonio and Collado-Vides, 2003; Perrenoud 
and Sauer, 2005). 

In addition to DNA-binding transcription factors, small 
regulatory RNAs, and other specific regulators, gene expres- 
sion also responds to changes in the overall physiological state 
of the cell. The latter comprises a variety of physiological 
parameters that have an impact on the rate of transcription and 
translation, such as the concentrations of (free) RNA 
polymerase and ribosome, gene copy number, and the size 
of amino acid and nucleotide pools. Contrary to specific 
regulators, these parameters affect the expression of all genes. 
Classical studies in bacterial physiology (Schaechter et aU 
1958; Kjeldgaard, 1961; Neidhardt and Fraenkel, 1961; Bremer 
and Dennis, 1996) reviewed by Maal0e and Kjeldgaard (1966) 
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and Scott and Hwa (2011), have demonstrated that in steady- 
state conditions the growth rate can be used to characterize 
this global physiological state. 

The joint control of gene expression by both specific effects 
of transcription factors and global effects of the physiological 
state has received relatively little attention thus far. Among the 
exceptions, we cite the work of Klumpp et al (2009), who have 
shown that the steady-state concentration of proteins in simple 
network architectures depends on the combined action of 
transcription factors and the growth rate. Dennis et al (2004) 
review the large amount of data on the control of rRNA 
synthesis in E. coli accumulated over several decades. They 
propose a model that integrates both growth rate-dependent 
effects and the regulatory control exerted by transcription 
factors. Notwithstanding the insights gained from these and 
other studies (Snoep et al 2002; Tadmor and Tlusty, 2008; Tan 
et aU 2009), they are limited in two respects. First, they 
consider the control of gene expression at steady state, not 
during transitions between physiological states. Second, there 
is currently no data set available that allows to study the 
contributions of transcription factors and global physiological 
effects at the level of a regulatory network. 

Here, we address the above questions in the context of 
carbon metabolism in Escherichia coli. We study a central 
regulatory circuit in this complex system, consisting of two 
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pleiotropic regulators of the cell: Crp (Gosset et aU 2004} and 
Fis (Bradley et aU 2007} . These transcription factors regulate, 
in response to the availability of carbon sources in the 
environment, a large number of genes encoding enzymes in 
central metabohsm (Gutierrez-Rios et aU 2007; Baldazzi et aU 
2010}. An example is the gene acs (Wolfe, 2005}. This gene 
encodes the enzyme acetyl-CoA synthetase (Acs}, which 
converts acetate to acetyl-CoA. acs is strongly expressed in 
the absence of glucose and is thus an excellent indicator of the 
transcriptional response of carbon metabohsm to a growth- 
phase transition. Figure 1 shows the dense regulatory structure 
of the network which, in addition to the transcription factors, 
involves cyclic AMP (cAMP} . This signaling metabolite forms 
a complex with Crp that binds to the promoter region of target 
genes (Kolb et aU 1993}. 

How does the regulatory circuit in Figure 1 respond to 
glucose depletion by dynamically adapting the expression of 
its genes? And what are the relative contributions of 
transcription factors and global physiological parameters 
to the observed changes in gene expression? To answer these 
questions in a quantitative way, we monitored the promoter 
activity of the genes in real time and in vivo, by means of GFP 
reporters. In parallel, a GFP reporter driven by a constitutive 
promoter was used to assay the time-varying physiological 
state. We show that a simple, parameterless mathematical 
model, in combination with careful data analysis procedures, 
can be used to separate the variation of the promoter activity of 
the genes into a part due to global physiological control and a 
part due to the effect of transcription factors. To verify if the 
latter part can be accounted for by known regulators, in 
particular Crp • cAMP, we extended the model and measured 
the time-varying concentration of cAMP. The above experi- 
ments were repeated under various physiological and genetic 
perturbations. 
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Figure 1 Central regulatory circuit involved in the control of £ coli carbon 
metabolism, consisting of the two pleiotropic transcription factors Crp and Fis. 
Crp is activated by the signaling metabolite cAMP, which accumulates upon the 
diminution of glycolytic fluxes. Crp • cAMP stimulates the expression of the gene 
acs, an effect counteracted by Fis (Wolfe, 2005). In addition, Crp • cAMP and Fis 
cross-regulate and auto-regulate the expression of their own genes as well as a 
large number of other genes (Hanamura and Aiba, 1 991 ; Ninnemann et al, 1 992; 
Ishizuka etal, 1994; Gonzalez-Gil etal, 1998; Nasser etal, 2001; Keseler etal, 
2011). While Fis inhibits the transcription of crp and its own gene, Crp • cAMP 
inhibits the transcription of fis and both activates and inhibits the transcription of 
crp, depending on its binding site in the promoter region. The global physiological 
state affects the expression of all genes in the network. Degradation and growth 
dilution balance protein synthesis. The effect of the growth rate on protein dilution 
has been omitted in order not to clutter the figure. Genes are shown in blue and 
promoters in red. Specific regulatory interactions are indicated by dashed lines 
and the effect of the physiological state by solid lines. 



The results of the above analysis provide novel insights into 
the response of the bacterial cell to nutritional stress. In 
particular, we show that the dynamic control of gene 
expression is shared between a signaling metabolite cAMP 
and the global physiological state of the cell, that is, the time- 
varying activity of the gene expression machinery and other 
global parameters. The DNA-binding transcription factors were 
not found to exert a significant control on gene expression. This 
observation is surprising, because Fis and Crp are located at the 
highest level of the hierarchical transcription regulation 
network in E. colt Their mutual regulatory interactions form 
a cross-inhibition motif, associated with bistability in synthetic 
circuits as well as in naturally evolved networks (Gardner et al, 
2000; Sanchez and Thieffry, 2001; Manu et al 2009; Graham 
et aU 2010} . It suggests a physiological role for the pattern of 
regulatory interactions, in that the circuit may function as a 
regulatory master switch controlling the adaptation of gene 
expression in response to carbon depletion and other stresses 
(Ropers et al 2006). 

Our data do not provide evidence that such a master switch 
is operative in the cell. More generally, they call for a 
reappraisal of the role of transcription factors in the coordina- 
tion of gene expression changes during growth transitions. 
Whereas transcriptional regulatory networks have sometimes 
been seen as the driving force of the genome-wide adaptation 
of gene expression, it may be more appropriate to view them as 
complementing the global control exerted by the physiological 
state of the cell. This shift in perspective, which places an 
often-neglected aspect of gene regulation in the foreground, 
has an obvious theoretical interest. However, it is also 
important in practice for the interpretation of transcriptome 
data (Regenberg et aU 2006} and the (re) design of biological 
systems in biotechnology and synthetic biology (Shachrai et aU 
2010; Scott et aU 2010; Carrera et al 2011). The method we 
present to quantify the relative contributions to the control of 
gene expression by specific and global physiological effects 
can be easily transposed to other regulatory systems in 
bacteria and higher organisms. 

Results 

Monitoring the dynamic response of the network 

To experimentally characterize the dynamic response of the 
network to glucose depletion, we systematically measured the 
input signals of the network in Figure 1, the global 
physiological state of the cell and the concentration of cAMP. 
In parallel, we monitored the outputs of the system, the time- 
varying activities of the acs, crp, and fis promoters (Figure 2} . 

Batch cultures of bacteria were grown in a microplate in 
minimal medium with glucose. We started the measurements 
after about 600 min (corresponding to about seven genera- 
tions) when the system had reached a well-defined initial state. 
The shape of the absorbance curves is typical for growth in 
minimal medium: exponential growth of the bacterial popula- 
tion, followed by a sharp drop of the growth rate due to glucose 
exhaustion (Figure 3 A) . 

At chosen time points along the growth curve, we 
determined the concentration of external cAMP using a 
luminescence-based immunoassay. From these 
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measurements, we estimated the internal cAMP concentration 
by means of a kinetic model accounting for cAMP import and 
export, as explained in Supplementary Section S3 and 
Supplementary Figures S3 and S4. The shape of the 
intracellular cAMP concentration profile agrees very well with 
other, direct measurements (Makman and Sutherland, 1965; 
Kao et aly 2004) . cAMP concentrations are low in the presence 
of glucose, rapidly accumulate at the end of exponential 
growth, when glucose is exhausted, and return to a lower 
steady-state level at the end of the transition (Figure 3B). 

In parallel, we monitored the promoter activities of the 
genes in the network using reporter plasmids carrying a 
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Figure 2 Real-time monitoring of gene expression. Bacteria recovered from 
glycerol stock are grown overnight in minimal medium with glucose, and then 
strongly diluted into the same medium in a 96-well microplate. When the culture 
has reached a quasi-steady state, after about 500 min, we monitor the growth of 
the bacterial population and the activity of a fluorescent reporter gene carried on 
a low-copy reporter plasmid. The primary absorbance and fluorescence data are 
treated by data analysis procedures, yielding the time-varying quantification of 
the growth rate and the promoter activities. Source data for this figure is available 
on the online supplementary information page. 



transcriptional fusion of a gfp reporter gene with the fis, crp, 
and acs promoters. The results are shown in Figure 4A-C. The 
promoter activities of fis and crp gradually decrease at the end 
of exponential phase, and then remain at a basal level after the 
exhaustion of glucose, with a slight recovery toward the end of 
the experiment. The latter may be a consequence of the 
consumption of acetate excreted during growth on glucose. 
The observed promoter activity of fis is in good agreement 
with northern blot quantifications of mRNA (de Jong et at, 
2010) . The activity pattern of acs in exponential phase seems to 
decrease as well, but as the fluorescence signal is close to the 
background level, the confidence intervals are wide and do not 
allow an unambiguous conclusion to be drawn. However, 
contrary to what was observed for fis and crp, the expression of 
acs is induced when glucose is exhausted at about 700 min. 
This observation is consistent with other reports in the 
literature (Wolfe, 2005). 

In steady-state conditions, the global physiological para- 
meters of the cell correlate with the growth rate (Bremer and 
Dennis, 1996} . The time-varying physiological state of the cell, 
such as the concentration of free RNA polymerase, is difficult 
to measure though (Klumpp and Hwa, 2008). As an indirect 
read-out of the global physiological state we therefore decided 
to use a constitutive promoter, whose activity is controlled by 
the transcription and translation machinery and the pools of 
precursor metabolites, but not by any particular transcription 
factor (Liang et aU 1999). 

As our constitutive promoter, we chose the pRM promoter of 
phage X. The transcription factors known to bind to this 
promoter, CI and Cro, are specific to the phage and are not 
present in uninfected E. coli cells (Oppenheim et aU 2005} . A 
priori, the activity of pRM therefore seems a good indicator of 
changes in the overall physiological state of the cell. As a 
control of the suitability of the choice of pRM, we repeated the 
experiment of Figure 4 with another promoter believed to be 
constitutively transcribed. In particular, we replaced pRM by 
the ptet promoter (Klumpp et aU 2009). The time-varying 
activities of the two promoters were observed to agree well 
in our conditions, both qualitatively and quantitatively 
(Supplementary Section S4; Supplementary Figure S5). This 
coincidence makes it unlikely that pRM is controlled by an 
unknown regulatory factor (as such a factor would have to 
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Figure 3 Experimental monitoring of physiological parameters. (A) Growth rate (• , blue) as computed from the measured absorbance of a bacterial culture (-, red). 
(B) Intracellular concentration of cAMP in wild-type strain ( • , blue) as derived from measured external concentrations of cAMP. (C) Idem for Ms strain. The data shown 
in the plots are the mean of 3-4 experimental replicates, with confidence intervals computed from the standard error of the mean (see Materials and methods). Source 
data for this figure is available on the online supplementary information page. 
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Figure 4 Experimental monitoring of transcriptional response of network. (A) Time-varying activity of fis promoter (#, blue), derived from GFP reporter data, and 
absorbance (solid line, red). (B-E) Idem for the activities of the crp, acs, and rpoS promoters, as well as the activity of the pRM promoter of phage X. The latter promoter is 
constitutive in our conditions and reflects the global physiological state of the cell. The primary fluorescence data for these curves are shown in Supplementary Section S9 and 
Supplementary Figure S9. (F) Idem for the activity of the acs promoter in a Afis strain. Source data for this figure is available on the online supplementary information page. 



impact both promoters in exactly the same way) . As a further 
control, we also excluded that cAMP has a significant 
regulatory effect on transcription from the pRM promoter 
(Supplementary Section S5; Supplementary Figure S6}. 

We thus conclude that variations in the activity of the pRM 
promoter reflect changes in the global physiological para- 
meters of the cell. This allows us to monitor the physiological 
state of the cell in real time and in vivo during the growth 
transition. The promoter activity computed from the fluores- 
cence signal is shown in Figure 4E. The activity of the pRM 
promoter is seen to be approximately stationary in exponential 
phase and stabilizing to a lower value after growth arrest. 

The advantage of plasmid-borne reporters is that they 
generally generate a strong signal, well beyond the autofluor- 
escence background. However, it is important to bear in mind 
that the plasmid copy number may vary with the growth rate 
(Lin-Chao and Bremer, 1986} and thus introduce a quantitative 
bias. We measured the variation of the number of plasmids per 
chromosomal equivalent of DNA across growth phases using 
quantitative PGR (qPCR), and found that it increases up to a 

4 Molecular Systems Biology 2013 



factor of 2 during the growth transition (Supplementary 
Section S6; Supplementary Figure S7) . This does not invalidate 
the qualitative shape of the profiles, especially the fall in 
activity of the constitutive promoter (which is actually 
underestimated) . However, to achieve quantitative precision, 
we need to develop an analysis method that corrects for this 
bias (see below). 

Dissection of the control of gene expression 

To analyze the relative contributions of transcription factors 
and the physiological state to the transcriptional response of 
the E. coli regulatory circuit, we use a simple mathematical 
model of promoter activity. Let pit) denote the promoter 
activity (M min~ ^) as a function of time t (min). We write 

P{t)=kp,{t)p2{t) (1) 

with k (Mmin~^) representing the maximum promoter 
activity. The dimensionless term pi(0, for convenience 
assumed to vary between 0 and 1, quantifies the modulation 

© 2013 EMBO and Macmillan Publishers Limited 



Shared control of gene expression in bacteria 
S Berthoumieux ef a/ 



of the promoter activity by the global physiological state, for 
instance through the availability of free RNA polymerase 
(Klumpp etaU 2009; Kotte etaU 2010). The dimensionless term 
Pzit)^ also varying between 0 and 1, accounts for the effect of 
transcription factors and other specific regulators, and may 
take the form of sigmoidal regulation functions often found in 
gene network modeling (de Jong, 2002; Bintu et aU 2005). 

To eliminate the unknown constant /c, we normalize 
Equation (1) with respect to a reference state at time In 
what follows, the reference state is chosen at the growth 
transition, indicated by the peak in acs activity, but other 
reference states are possible as well. We define p^=p(f), 
pO=pi(tO), and p\=P2{f), from which it follows with 
Equation (1) thatp° = kp\p2. Division oip(f) by p° gives, after 
a logarithmic transformation. 



lO! 



Pi 



(2) 



Two special cases of this model, which allow to answer two 
questions of specific biological interest, will be examined in 
more detail. 

(i) When the global physiological effect is dominant, that is, 
when the effect of the transcription factors is negligible, we 
have p2{t) ^ P2 and the second term in the right-hand side of 
Equation (2) approximates 0. Bearing in mind that the global 
effect is measured by the activity of the constitutive pRM 
promoter, we can rewrite the model as 



(3) 



with Prm(0 and p^j^ the time-varying activity of the pRM 
promoter and its value at respectively. 

(ii) The data in Figures 3 and 4 indicate that both the 
promoter activity of acs and the intracellular concentration of 
cAMP peak after growth arrest. This suggests a simphfied 
model for acs, and potentially other Crp • cAMP-regulated 
genes, in which the remaining variation of the promoter 
activity after subtraction of global physiological effects is 
proportional to the intracellular cAMP concentration c{t) [M] : 



y Prm ^ 



(4) 



with c° = c{t^) . Notice that this model is based on the data, but 
that the biological assumptions underlying the simphfication 
can be exphcitly formulated (Supplementary Section S7). 

As will be shown below, the models of Equations (3) and (4) 
can be straightforwardly tested in different conditions by 
means of experimental data, by inserting for p(0 the measured 
promoter activities of fis, crp, and acs (denoted by Pfisit], 
Pcrp(0, andpacs(0» respectively). This allows us to answer the 
following questions, (i) To which extent can the observed 
variation in the promoter activity of the genes in the network of 
Figure 1 be accounted for by the effect of the global 
physiological state only? (ii) How much of the remaining 
variation in promoter activity can be explained by changes in 
the intracellular concentration of cAMP? 

Notice that the models of Equations (3) and (4) have a 
number of advantages for this purpose. First, they are 
parameterless and therefore do not require preliminary model 
calibration. Second, when the same plasmid vector is used for 



measuring the activity of a constitutive and a network gene, 
the growth phase-dependent variation of the plasmid copy 
number equally affects the target promoter p{t) and the 
constitutive promoter Prm(0, and therefore cancels out in the 
equations. If different vectors are used, then the equations can 
be easily adapted to correct for this bias by means of qPCR data 
(see Supplementary Section S8 and Supplementary Figure S8 
for the mathematical arguments and an illustration) . 

Shared control of gene expression by 
transcription factors and the global physiological 
state of the cell 

We first tested the hypothesis that the adaptation of gene 
expression in response to glucose exhaustion is mainly 
controlled by the physiological state of the cell, as measured 
by the activity of the pRM promoter. In this case. Equation (3) 
predicts a linear relation between log{p{t)/p^) and 
log(pRM(0/PRM)' diagonal in the scatter plots of 

Figure 5A-C. If the global effects are dominant, then one 
would expect the data points to be scattered around the 
diagonal. As can be seen, this is indeed the case for fis and crp, 
but not for acs. 

To quantify the fit of the models with the data, we computed 
the coefficient of determination (i^^), the square of the 
correlation coefficient (Hamilton, 1992). For the computation 
of the correlation coefficient, the data points were weighed by 
the square inverse of the size of the confidence interval, 
ensuring that data points with low uncertainty weigh more 
than points with high uncertainty. The coefficient of determi- 
nation measures the proportion of the variance predicted by the 
model. For fis and crp, we have high values (0.93 and 0.96). 
However, for acs the value is found to be much lower (0.08) . 

To account for the unexplained variation of the promoter 
activity of acs, we extended the model with the known 
regulation of this gene by Crp • cAMP. Equation (4) predicts a 
linear relation between the variation of acs activity remaining 
after subtraction of the global effect, and the intracellular 
concentration of cAMP. This prediction corresponds to the 
diagonal in Figure 5D. The experimental data points are seen 
to correspond well with the predictions of the extended model 
(i^^ = 0.93). We also tested if the addition of Crp • cAMP 
regulation of fis and crp could account for the remaining 
variance in the data. This turned out not be the case, in the 
sense that we obtained very low R^ values when extending the 
model with an additional regulator (0.03 and 0.02, respec- 
tively, see Figure 51 for an example). We therefore conclude 
that the time-varying variation of fis and crp activity is well 
accounted for by changes in the physiological state, whereas 
for acs we also need to include the accumulation of cAMP. 

Another informative view on the same data can be obtained 
by making explicit the contribution of specific regulatory 
factors to the observed variation in promoter activity. Equation 
(1) impHes that 



p(0 ^ PrmW Pijt) 
p'km p\ ■ 



(5) 



As a consequence, the (normalized) contribution of specific 
regulatory factors can be obtained by dividing the (normalized) 
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Figure 5 Predicted and observed control of fis, crp, and acs activity by Crp • cAIVIP and the physiological state of the cell, in various experimental conditions and genetic 
backgrounds. (A) Predicted (-, black) and measured (#, blue) relative activity of the fis promoter (log(p^is(t)/p°^J) as a function of the relative activity of the pRM 
promoter (log(pRM(0/PRM))- The confidence intervals in the plots have been computed from experimental replicas, as described in Materials and methods. 
(B, C) Idem for crp and acs. (D) Predicted (-, black) and measured (•, blue) remaining relative activity of the acs promoter after subtraction of the effect of global 
physiological parameters (lo g(pacs(0/Pacs) ~ 1o8(Prm(0/Prm)) ^"^^ ^ function of the relative intracellular cAMP concentration (log(c(t)/c°)). (E) Same as (A), 
but in a Acrp strain. (F-H) Same as (B-D), but in a Afis strain. (I) Same as (D), but for crp. (J-L) Predicted (-, black) and measured ( • , blue) relative activity of the rpoS 
promoter (log (prpos(0/P?pos)) ^ function of the relative activity of the pRM promoter (log(pRM(0/PRM)) '"^ ^^^^^ different conditions (wild-type, Afis, and Acrp). 
Source data for this figure is available on the online supplementary information page. 



promoter activity of the gene of interest by the (normalized) 
promoter activity of the constitutive gene. Figure 6 shows the 
variation of the ratio ^2(^/^2 ^^^^ duration of the 
experiment for the genes in our network. As expected, the 
ratio remains close to 1 for fis and crp, while for acs it follows 
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the intracellular concentration of cAMP (Figure 3B} . Some small 
deviations of this pattern can be seen though. For example, while 
overall the acs promoter behaves as a constitutive promoter, 
there seems to remain a small, regulated component in its 
dynamics before the growth transition. This conclusion should 
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Figure 6 Relative contributions of specific regulators to the time-varying activity of the network genes. (A) Relative promoter activity p^i5(f)/p^.^ (-, blue) and relative 
effect of transcription factors and other specific regulators p^js,2 (0/Pfis 2 ("> P''^'^)- ^) ^^"^^ ^''P ^"^^ Pfis,2 Wp% 2 Pcrp,2 (0/P?rp,2 ■'©'^ain 

close to 1 , indicating that the effect of specific regulation is negligible for ihese promoters. The variation of pacs,2 (0/Pacs,2 corresponds to the cAMP concentration profile 
(Figure 3B), suggesting a regulatory role for this signaling metabolite. Source data for this figure is available on the online supplementary information page. 



be taken with care though, as the error bars are large due to 
difficulties in treating the fluorescence signal emitted by the acs 
reporter (which is close to the autofluorescence background 
during exponential growth). 



Validation of predicted dominance of global 
physiological effects 

The surprising observation that growth phase-dependent 
effects dominate the expression dynamics of fis and crp calls 
for further experimental tests. To verify if changes in the 
medium after extended growth on glucose might account for 
the results, we performed a down-shift experiment. Once the 
culture reached quasi-steady-state growth in minimal medium 
with 0.3% glucose, we rediluted the bacteria into the same 
medium with a low glucose level (0.06%), thus enforcing a 
rapid growth arrest. In parallel, we monitored the expression 
of the genes in the network and the activity of the pRM 
promoter. The results are shown in Supplementary Sections S9 
and SIO; Supplementary Figures S9-S14; Supplementary Table 
S4. They confirm, both qualitatively and quantitatively, the 
observations made in the reference conditions: the variation in 
the activity of the fis and crp promoters is mainly controlled by 
the global physiological state, whereas the variation in the 
promoter activity of acs is to a large extent predicted by the 
combined effect of the growth phase and cAMP. 

If the control exerted by the physiological state of the cell 
accounts for the major part of the variation of fis and crp 
activity, that is, if the cross-regulatory and autoregulatory 
interactions have a minor effect in our conditions, then one 
would expect the minimal model of Equation (3) to predict the 
variation in the promoter activity equally well in Afis and Acrp 
backgrounds. To test this prediction, we measured the input- 
output behavior of the network in mutant strains deleted for fis 
and crp. The resulting data were used to construct Figure 5E 
and F. The plots confirm the prediction that the global 
physiological effect is dominant in the control of the expression 
of crp in a Afis strain and fis in a Acrp strain {R^ = 0.70 and 
R^ = 0.96, respectively). That is, cross -regulation between Fis 
and Crp is of little importance here. This also holds for 
autoregulation, since the dynamics of crp activity in a Acrp 
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strain and fis in a Afis strain are also well accounted for by the 
model of Equation (3) [R^ = 0.98 and R^ = 0.95, respectively) . 

Similarly, in the case of acs, one would expect the results of 
the analysis to change little in a Afis mutant, given that the 
overall cellular physiology and cAMP were found to dominate 
its expression control. This prediction is also confirmed by our 
data, as shown in Figure 5G and H. The effects of Crp • cAMP 
and the global physiological state together account for most of 
the variation in acs activity (i^^ = 0.96), whereas global effects 
alone fail to be a good predictor [R^ = 0.56) . On the other hand, 
the deletion of crp disables the sensor mechanism of glycolytic 
fluxes provided by Crp • cAMP (Bettenbrock et at, 2007), and 
is known to prevent the induction of acs when glucose is 
exhausted (Wolfe, 2005) . This is confirmed in our data: acs is 
not expressed in a Acrp background (Supplementary Section 
S9; Supplementary Figure Sll). 

The results obtained thus confirm the observation that 
growth phase-dependent effects dominate the transcriptional 
control of the pleiotropic regulators Fis and Crp. The question 
can be raised to which extent this holds for other key 
transcriptional regulators involved in the response to glucose 
depletion. We addressed this question for RpoS or a^, a master 
stress regulator of E. coli and other bacteria (Hengge-Aronis, 
2002; Battesti et aU 2010) . The main rpoS promoter is contained 
within the nlpD gene and is negatively regulated by Crp • cAMP 
(Lange and Hengge-Aronis, 1994). The regulatory mechanism 
is not well understood and the effect remains somewhat 
controversial (see Zgurskaya et al, 1997 and references 
therein) . In the same way as for the other genes, we monitored 
the expression of a transcriptional fusion of the rpoS promoter 
with a gfp reporter gene (Supplementary Section S9; 
Supplementary Figure SI 3). The data thus obtained were 
analyzed by means of the models of Equations (3) and (4), 
giving the results shown in Figure 5J-L. As for fis and crp, the 
transcription of rpoS follows the activity of the constitutive 
pRM promoter quite well, in all considered conditions. 

Discussion 

The variation of gene expression across growth phases is 
controlled both by the physiological state of the cell and by 
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transcription factors. Here, we have shown how to distinguish 
between these two effects, using a simple mathematical model 
of promoter activity and carefully designed data analysis 
procedures. The approach has several advantages that make it 
easy to put to work in bacteria but also in higher organisms. 
The models do not have free parameters that need to be 
calibrated, hypotheses on the effect of regulators can be readily 
tested by monitoring the expression of target genes and a 
constitutive control, and the use of plasmid-borne reporter 
systems does not bias the analysis. The main novelty is that 
this allows the relative contributions of specific regulators and 
the global state of the cell to be monitored in vivo, both 
dynamically in time and on the level of a regulatory network. 

We apphed the method to expose hidden correlations in a 
circuit involving key regulators of carbon metabolism in 
E. colt A first surprising finding of this study is that the activities 
of the genes encoding the transcription factors Fis, Crp, and 
RpoS closely follow the activity of a constitutively transcribed 
promoter. The activity of this promoter is assumed to reflect the 
global physiological state of the cell, that is, parameters affecting 
the expression of all genes, such as the concentrations of (free) 
polymerase and ribosome, gene copy number, and the size of 
the pools of amino acids and nucleotides. 

An alternative interpretation of the observed correlations 
would be that a specific stress factor be responsible for the 
simultaneous adjustment of all above-mentioned promoter 
activities. For example, guanosine 3',5'-bispyrophosphate or 
(p)ppGpp (Potrykus and Cashel, 2008) could have this role. 
This interpretation is motivated by the fact that the transcrip- 
tion of some of the genes in the network of Figure 1 is under 
stringent control and the fis promoter is directly regulated by 
(p)ppGpp (Ninnemann et al, 1992; Walker et al, 1999; 
Johansson et al, 2000; Mallik et al, 2004). For this explanation 
to hold, however, we have to postulate that (p)ppGpp has 
exactly the same specific, quantitative effect on all genes in 
our network as well as on the two constitutive promoters 
tested (pRM and ptet, see Supplementary Section S4 and 
Supplementary Figure S5). This is unlikely to be the case, so 
we attribute the observed variations in promoter activity to 
changes in the global physiological state of the cell. Notice that 
the absence of a (strong) specific effect of (p)ppGpp, at least in 
our conditions, does not contradict stringent control of the 
network genes. (p)ppGpp is a major factor in the control of the 
global physiological state (Traxler et al, 2006; Potrykus et al, 
2011). It inhibits transcription of the rRNA operons, activates 
amino-acid biosynthesis promoters, and indirectly influences 
the availability of free RNA polymerase by inhibiting strong 
promoters (Durfee etal, 2008; Potrykus and Cashel, 2008), 
thus influencing the activity of the gene expression machinery. 
Through these mechanisms, (p)ppGpp may have an indirect 
effect on the expression of a large number of genes. In this 
study, we have captured these indirect, global effects into an 
easy-to-measure variable, the activity of a constitutive gene. 
An interesting extension of the model would be to analyze the 
global control of promoter activities in more detail, distin- 
guishing between the contributions from individual physiolo- 
gical parameters. 

A second interesting observation, actually a consequence of 
the first, is that in our conditions the regulatory interactions 
involving Fis and Crp do not significantly contribute to the 
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Figure 7 Reduced regulatory circuit. Tlie network includes the interactions that 
were found to dominate the transcriptional response of the network in Figure 1 , in 
particular the activation of all genes by the physiological state of the cell and the 
activation of acs by Crp • cAMP. 



control of the network response (Figure 7). Fis and Crp are 
among the most pleiotropic transcription factors of the cell, 
located at the top of the hierarchically structured transcription 
regulation network of E. coli (Martinez-Antonio and Collado- 
Vides, 2003). Their regulatory interactions form a cross- 
inhibition motif, which suggests that the network behaves as a 
bistable switch, along the lines of what has been shown for 
synthetic circuits in bacteria and gene networks underlying 
developmental processes in higher organisms (Gardner et al, 
2000; Sanchez and Thieffry, 2001; Manu et al, 2009; Graham 
et al, 2010) . It notably raises the question if the adaptation of 
gene expression in response to glucose depletion might be 
coordinated at the highest level of the transcription regulation 
hierarchy by a master switch under the control of cAMP 
(Ropers et al, 2006). Our data do not support this transcription 
factor-oriented view of network functioning, however, in the 
sense that the observed expression dynamics are very well 
explained without taking into account the dense pattern of 
regulatory interactions between Fis and Crp. Rather, the 
driving force behind the adaptation of gene expression is a 
global physiological effect that usually does not even figure in 
regulatory network diagrams. 

The absence of clearly identifiable effects of Fis and 
Crp • cAMP, other than on acs, is surprising as these interac- 
tions have been well documented in the literature (Ishizuka 
etal, 1994; Gonzalez-Gil etal, 1998; Nasser etal, 2001) and are 
beheved to be important in the adaptation to glucose 
exhaustion. Moreover, Fis and Crp have a large number of 
evolutionary-conserved binding motifs in the E. coli genome 
(Grainger et al, 2005; Cho et al, 2008). This paradox may be 
resolved by observing that when the physiological role of the 
above-mentioned interactions was tested, for example the 
activity of the fis promoter in a Acrp strain (Nasser et al, 2001) 
or the activity of rpoS in the same background (Zgurskaya et al, 
1997), no dramatic effects were found. Rather than dominating 
the control of gene expression, the specific regulators seem 
to finetune a global trend set by the availability of the 
transcriptional and translational machinery and other global 
physiological parameters. This finetuning may nevertheless be 
important for the cell, as even small regulatory effects may 
confer a decisive growth advantage. For instance, experimen- 
tal evolution studies in E. coli have shown that mutations in a 
transcription factor providing a fitness benefit of <1% 
outcompete isogenic strains not carrying this mutation 
(Pelosi et al, 2006). 
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While the control of gene expression during growth 
transitions is thus shared between global physiological effects 
and specific transcription factors, our results question the 
central role often attributed to transcriptional regulatory 
networks in controlling genome-wide expression changes 
during physiological transitions. It may be more appropriate 
to regard transcriptional regulators as complementing and 
finetuning the global control exerted by the physiological state 
of the cell. It is therefore strongly recommended to exphcitly 
take into account these global parameters when analyzing 
transcriptome data or when designing a synthetic controller 
circuit. 

Materials and methods 
strains and growth conditions 

The E. coli strains used in this study are the wild-type strain BW25113 
and the deletion mutants Afis and Acrp. The wild-type and mutant 
strains were transformed with low-copy plasmids bearing a fusion of a 
gfp reporter gene with the promoter regions of the genes crp, fis, acs, 
and rpoS (de Jong et al, 2010). In addition, we used a plasmid carrying 
the pRM promoter of phage X and the synthetic ptet promoter. More 
details on the strains and plasmids used in this study can be found in 
Supplementary Section SI and Supplementary Tables SI -S3 of 
Supplementary information. 

The strains were recovered from glycerol stock and grown overnight 
at 37° C in M9 minimal medium (Miller, 1972) supplemented with 
0.3 % glucose and mineral trace elements. The overnight cultures were 
strongly diluted (1500- to 7000-fold) into a 96-well microplate, so as to 
obtain an adjusted initial ODeoo of 0.001. The wells of the microplate 
contain 150 jil of the above medium, to which was added 1.2% of the 
buffering agent HEPES (4-(2-hydroxyethyl)-l-piperazineethanesulfo- 
nic acid) for maintaining a constant external pH. The wells were 
covered with 60 |il of mineral oil to avoid evaporation. The microplate 
cultures were then grown for 20h at 37°C, with agitation at regular 
intervals, in a microplate reader (Fusion Alpha, Perkin-Elmer) . 

Real-time monitoring of gene expression and data 
analysis 

The expression of the reporter genes in the different genetic 
backgrounds was monitored in real time. During a typical acquisition 
period, we obtain about 120 readings each of absorbance (600 nm) and 
fluorescence (485-520 nm). The promoter activities were computed 
from these data by means of carefully designed data analysis 
procedures implemented in Matlab. Extending earlier work (Ronen 
et al, 2002; de Jong et al, 2010), we notably took into account the effect 
of protein degradation, maturation time, dynamic changes of 
autofluorescence, and synchronization of experimental replicates 
(see Supplementary Section S2 and Supplementary Figures SI and 
S2 for details) . Confidence intervals consist of ± 2 standard errors of 
the mean, determined from on average four experimental replicas. 
Under the assumption of Gaussian distributions, this corresponds to 
95% confidence intervals. 



Measurement of cAMP concentrations 

To measure the extracellular concentration of cAMP (adenosine 3',5'- 
cyclic monophosphate), we used a commercially available immu- 
noassay kit (Upstate). We took 100 |il samples at regular time intervals 
from cultures growing in a microplate, under the conditions described 
above (12 time points, 3 replicates). Following the manufacturer's 
instructions, the cAMP concentration at the different time points was 
determined from luminescence measurements in the microplate 
reader and a calibration standard relating luminescence intensity to 
cAMP concentration. 



Measurement of time-varying plasmid copy 
number 

We used qPCR to determine the time-varying number of plasmids per 
chromosomal equivalent of DNA (plasmid copy number), following a 
previously validated protocol (Lee et al, 2004). Details of the 
experimental procedure can be found in Supplementary Section S6 
and Supplementary Figure S7. 



Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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